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I. INTRODUCTION 



In this paper we study the process 

e + e~ e + e~7/ c , (1) 
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where the r\ c meson (J PC = h ), the lowest lying char- 
monium state, is produced via the two-photon produc- 
tion mechanism illustrated in Fig. [TJ We measure the 
differential cross section for this process in the single-tag 
mode where one of the outgoing electrons 1 is detected 
(tagged), while the other electron is scattered at a small 
angle and hence is undetected (untagged). The tagged 
electron emits a highly off-shell photon with squared mo- 
mentum transfer q 2 = —Q 2 = (p' — p)' 2 , where p and 
p' are the four-momenta of the initial- and final-state 
electrons. The momentum transfer squared to the un- 
tagged electron (gf) is near zero. The differential cross 
section da(e + e~ — > e + e~P)/dQ 2 for pseudoscalar meson 
(P) production depends on only one form factor F(Q 2 ), 



Unless otherwise specified, we use the term "electron" for cither 
an electron or a positron. 
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FIG. 1: The Feynman diagram for the e + e e + e rj c two- 
photon production process. 



which describes the 77* — > P transition. To relate the 
differential cross section to the transition form factor we 
use the formulae for the e + e~ — > e + e~ir° cross section in 
Eqs. (2.1) and (4.5) of Ref. [J. 

According to perturbative QCD (pQCD), the transi- 
tion form factor can be presented as a convolution of a 
calculable hard scattering amplitude for 77* — > cc with a 
nonperturbative light-cone wave function (/), )e 0. The 
measurement of the form factor allows us to test the 
pQCD prediction and to obtain information on the shape 
of the r\ c wave function. The Q 2 dependence of the form 
factor is studied theoretically in Refs. 0,13] using pQCD, 
and in Ref. [j| using the lattice QCD approach. 

The T] c transition form factor was measured by the L3 
Collaboration Q with very small data sample. In this 
paper we study the e + e~ — > e + e~r\ c reaction in the Q 2 
range from 2 to 50 GeV 2 . The rj c is observed via the r\ c — > 
KsK ± tt t decay, 2 which allows the r\ c to be selected with 
relatively low background. The rj c two-photon width and 
branching fractions are not well measured Q. Therefore, 
we also study no-tag data (Q 2 ~ 0), measure the product 
r(^c 77)S(?/ C — > KKtt), and normalize the transition 
form factor F(Q 2 ) to -F(O). The measured values of the 
rj c mass and width obtained in different experiments have 
a large spread ■ The high statistics no-tag data sample 
allows us to extract precise values for these parameters. 



II. THE BABAR DETECTOR AND DATA 
SAMPLES 

We analyze a data sample corresponding to an inte- 
grated luminosity of 469 fb _1 recorded with the BABAR 
detector Q at the PEP-II asymmetric-energy storage 
rings at SLAC. At PEP-II, 9 GeV electrons collide with 
3.1 GeV positrons to yield an e + e~ center-of-mass (cm.) 
energy of 10.58 GeV (the T(4S) resonance). About 10% 



of the data used in the present analysis were recorded 
about 40 MeV below the resonance. 

Chargcd-particlc tracking is provided by a five-layer 
silicon vertex tracker (SVT) and a 40-layer drift cham- 
ber (DCH) operating in a 1.5-T axial magnetic field. The 
transverse momentum resolution is 0.47% at 1 GeV/c. 
Energies of photons and electrons are measured with a 
CsI(Tl) electromagnetic calorimeter (EMC) with a reso- 
lution of 3% at 1 GeV. Charged-particle identification is 
provided by specific ionization measurements in the SVT 
and DCH, and by an internally reflecting ring-imaging 
Chcrcnkov detector (DIRC). Electron identification also 
uses shower shape in the EMC and the ratio of shower 
energy to track momentum. Muons are identified in the 
solenoid's instrumented flux return, which consist of iron 
plates interleaved with resistive plate chambers. 

The signal process is simulated with the Monte Carlo 
(MC) event generator GGResRc. It uses the formula 
for the differential cross section for pseudoscalar meson 
production from Ref. [l|. Two samples of signal events 
are produced: one for no-tag measurement without any 
kinematic restrictions, and the other with the restrictions 
on the momentum transfer values to the electrons Q = 
~q\ > 1.5 GeV 2 and -q 2 < 1 GeV 2 . The restriction on 
Q 2 for the tagged electron corresponds to the detector 
acceptance. The experimental criteria providing these 
restrictions for data events will be described in Sec. IIIII 
In the simulation of no-tag events we use the form factor 

F(a 2 a 2 \ m ° ] (2) 

expected in the vector dominance model. The form factor 
is fixed to the constant value F(0, 0) in the simulation 
of single-tag events. The produced r\ c decays into the 
KsK^nT final state. The simulation uses a three-body 
phase space distribution to describe this decay. 

The GGResRc event generator includes next-to- 
lcading-order radiative corrections to the Born cross sec- 
tion calculated according to Ref. [§]. In particular, it 
generates extra photons emitted by the initial- and final- 
state electrons. The formulae from Ref. @ are modi- 
fied to take into account the hadron contribution to the 
vacuum polarization diagrams. The maximum energy of 
the photon emitted from the initial state is restricted by 
the requirement 3 E* < 0.05i/s, where y/s is the e + e~ 
cm. energy. The generated events are subjected to de- 
tailed detector simulation based on GEANT4 [TcJ, and 
are reconstructed with the software chain used for the 
experimental data. Temporal variations in the detector 
performance and beam background conditions arc taken 
into account. 



3 Throughout this paper an asterisk superscript denotes quantities 
2 The use of charge conjugate reactions is implied throughout un- in the e+e - cm. frame. In this frame the positive z-axis is 

less explicitly stated otherwise. defined to coincide with the e — beam direction. 



FIG. 2: The Feynman diagrams for the background processes (a) e + e — > J/V>7, and (b) e + e — > e + e J/tp. 



The processes with a J/ip in the hnal state (Fig. [2]) , 
with J /if) decaying into 7^7, can imitate the process 
under study. The initial state radiation (ISR) process 
(Fig. [U(a)) contributes to the no-tag mode, while the 
J/ip bremsstrahlung process (Fig.[2Jb)) contributes back- 
ground to the single-tag mode. We simulate both pro- 
cesses with J/tp decaying to KsK ± tt^ z and also to i] c -f 
followed by r\ c — > KsK^n 11 . To estimate a possible back- 
ground from other two-photon processes we simulate the 
reaction e + e _ — > e + e~77 c 7r°. These events are generated 
with an isotropic distribution of the final state mesons in 
the rj c ir° rest frame. 



III. EVENT SELECTION 

We select e + e~ — >• e + e~r/ c candidates in the no-tag and 
single-tag modes, with zero and one detected electron, 
respectively. The decay mode r] c — > KsK^-k^ (K s — >• 
7r + 7r~) is used to reconstruct r/ c meson candidates. 

Events are selected with at least four (five for single- 
tag mode) charged-particle tracks. Since a significant 
fraction of events contain beam-generated spurious track 
and photon candidates, one extra track and any num- 
ber of extra photons are allowed in an event. The tracks 
corresponding to the charged kaon and pion must be op- 
positely charged, and must extrapolate to the interaction 
region. The kaon is required to be positively identified, 
while the pion track is required to be inconsistent with 
the kaon hypothesis. 

The track identified as an electron must originate from 
the interaction region and be in the polar angle range 
0.387 < 9 < 2.400 in the lab frame (0.64-2.69 in the 
e + e~ cm. frame). The latter requirement is needed for 
good electron identification. To recover electron energy 
loss due to bremsstrahlung, both internal and in the de- 
tector material before the DCH, the energy of any EMC 
shower close to the electron direction is combined with 
the measured energy of the electron track. The resulting 
cm. energy of the electron candidate must be greater 
than 1 GeV. 

A Ks candidate is formed from a pair of oppositely 
charged tracks fitted to a common vertex, and yield- 



ing an invariant mass value in the range 487.5-507.5 
MeV/c 2 , when the charged-pion mass is assigned to each 
track. The candidate is then refitted with a Ks mass con- 
straint to improve the precision of its momentum mea- 
surement. To suppress combinatorial background, the 
angle between the Ks candidate momentum and the line 
connecting its production and decay vertices (ipK s ) is re ~ 
quired to satisfy cos-0x s > 0.95. 

An lie candidate is formed from Ks, K^, and 
candidates fitted to a common vertex. The prelimi- 
nary selection criterion for no-tag events requires that 
I cos#*J > 0.95, where 0* o is the polar angle of the can- 
didate rj c in the e + e~ cm. frame. Figure [3] shows the dis- 
tribution of the invariant mass of the pions forming a Ks 
candidate {M^) for events satisfying this criterion. The 
shaded histogram demonstrates the effect of the require- 
ment cosipKs > 0.95. The transverse momentum of the 
r) c candidate in the e + e~ cm. frame is restricted to the 
range p*^_ < 0.25 GcV/c. The invariant mass distribution 
for r) c candidates is shown in Fig. 01 where for events with 
more than one r\ c candidate (about 0.4% of signal events), 
the candidate with the smallest value of p* ± is selected. 
The rjc peak from two-photon production and the J/ip 
peak from the ISR process e + e~ — > J/ipj are clearly seen 
in the invariant mass distribution. The shaded histogram 
shows the distribution for candidates rejected by the re- 
quirement p* ± < 0.25 GeV/c. This requirement limits the 
momentum transfers to the electrons. The value of the 
effective threshold for qf (i = 1, 2) is determined from the 
dependence of the detection efficiency on max(— qf, —q\) 
and is about 0.1 GeV 2 . Such a low q 2 threshold yields a 
model-independent extraction of F(0,0) from the no-tag 
data. We note that d<r/dq 2 dq 2 cx l/{q 2 q 2 ) at small \qf\ 
and \q 2 \. 

For single-tag events we combine an rj c candidate with 
an electron candidate and require | cos 8l Vc \ > 0.95, where 
0* rja is the polar angle of the momentum vector of the 
eri c system in the e + e~ cm. frame. The transverse mo- 
mentum of the erj c system is restricted to p* ± < 0.25 
GeV /c. The p\ distribution for data candidates is shown 
in Fig. [51 where the shaded histogram is the correspond- 
ing distribution for simulated signal events. The condi- 
tion on p*^_ limits the value of the momentum transfer 
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FIG. 3: The Mi^ distribution for Ks candidates in the no- 
tag data sample. The shaded histogram shows events rejected 
by the requirement cos ipK s > 0.95. The arrows indicate the 
region used to select event candidates. 
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FIG. 5: The transverse momentum distribution for erj c data 
candidates. The shaded histogram is for simulated signal 
events. Data candidates for which p* ± < 0.25 GeV/c (indi- 
cated by the arrow) are retained. 
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FIG. 4: The distribution of the invariant mass of r] c candi- 
dates in the no-tag data sample. The shaded histogram shows 
events rejected by the requirement p* ± < 0.25 GeV/c. 



to the untagged electron The effective q 2 threshold 
determined from simulation is about 0.1 GcV 2 . 

The emission of extra photons from the electrons in- 
volved leads to a difference between the measured and 
actual values of Q 2 . In the case of ISR, Q 2 



) 2 

• mcas 



( 3true(l + r 7)j where r 7 = 2E*/y/s. To restrict the energy 
of the ISR photon we use the parameter 



r s — E, 



\P*er, 



(3) 



where E* er] ^ and p* er)c are the cm. energy and momentum 
of the detected er\ c system. In the ISR case this parame- 
ter coincides with r 7 defined above. The r distributions 
for data and simulated signal are shown in Fig. [6] Candi- 
dates with —0.02 < r < 0.03 are retained. We note that 
the condition on r ensures compliance with the restric- 
tion r-y < 0.1 used in the MC simulation. 



FIG. 6: The distribution of r defined in Eq.© for data 
(points with error bars) and simulated signal events (shaded 
histogram). The arrows indicate the region used to select 
candidate events (-0.02 < r < 0.03). 



For two-photon events with a tagged positron (elec- 
tron), the momentum of the detected ery c system in the 
e + e~ cm. frame has a negative (positive) z-componcnt, 
while events resulting from e + e~ annihilation are pro- 
duced symmetrically To suppress the e + e~ annihilation 
background, event candidates with the wrong sign of the 
momentum z-component are removed. 

The distribution of the invariant mass of r\ c candidates 
for single-tag events satisfying the selection criteria de- 
scribed above is shown in Fig. [7] For events with more 
than one rj c candidate, the candidate with smallest p*^ is 
selected. Signals corresponding to r\ c and J/ ip production 
are seen clearly in the mass spectrum. The J/tp events 
are from the process e + e~ — > e + e~J/ip (see Fig. HJb)). 
The shaded histogram in Fig. [7] shows the distribution 
for candidates with the wrong sign of the er\ c momentum 
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FIG. 7: The invariant mass distribution for single-tag r} c can- 
didates. The shaded histogram is for events with the wrong 
sign of the erj c momentum z-component. 



FIG. 8: The KsK ± -k^ 1 invariant mass spectrum for no-tag 
data events satisfying (dashed histogram) and not satisfy- 
ing (solid histogram) the condition for ISR event selection 
(Eq.flU). 



z-componcnt. Since the numbers of events from e + e~ 
annihilation in the wrong- and right-sign data samples 
are expected to be approximately the same, this shows 
that the background from e + e~ annihilation peaking at 
the rj c mass is small. 



IV. 



FITTING THE KsK^ MASS SPECTRUM 
FOR NO-TAG EVENTS 



The mass spectrum for no-tag events exhibits the r\ c 
and J /if) peaks corresponding to the two-photon and ISR 
production mechanisms, respectively. The cm. momen- 
tum p* of the KsK^tt^ system for ISR events is equal 
to (^/i/2)(l-Mls. /s). In the mass region under study 
the detector acceptance strongly limits the efficiency for 
ISR events. Due to the asymmetry of the acceptance 
most of the detected ISR events have positive cos 6** c . It 
follows that the ISR events can be selected by requiring: 

p*/(l - M 2 KR Js) > 5.1 GeV/c, cos 6* Vc < 0. (4) 

The mass distribution for the events satisfying this condi- 
tion is shown by dashed histogram in Fig. [8] The selected 
event sample contains mostly ISR events with very little 
two-photon rj c admixture. The KsK^it^ mass distribu- 
tion for the events not satisfying Eq.(U]) is shown by the 
solid histogram in Fig. [8j The remaining J /if) events are 
from the ISR process with more than one photon emitted 
from the initial state. 

To determine the r/ c mass and width, and the number 
of events containing an 7y c , a binned likelihood fit is made 
to the distributions in Fig. [8] using a function consisting 
of signal (rj c and J /if) and background distributions. The 
bin size used in the fit is chosen to be 2.5 MeV/c 2 . The 
J /if) line shape is represented by the detector resolution 
function for ISR events. The r\ c line shape is described 
by a Breit-Wigner function convolved with the detector 



resolution function corresponding to two-photon produc- 
tion. In each case, the detector resolution function is 
obtained using MC simulation of the detector response. 
We use the nonrclativistic Breit-Wigner form 



(r/2) 2 



(M - Mf + (r/2) 2 



(5) 



where M is the KsK ± t: t invariant mass, and M and T 
are the r\ c mass and width. The changes in the values of 
the parameters, if a relativistic Breit-Wigner function is 
used, are negligible. 

The detector resolution functions are determined from 
the distributions of the difference between measured and 
true simulated KsK ±r K^ mass for the processes e + e~ — s- 
e + e~?/ c , 7] c — > I<sK ± it t and e + e~ — > J/ipj, J /if) — > 
KsK ± n T shown in Fig.^a) and Fig. [§Ib), respectively. 
The distributions are fit with the following function: 



where 



F{x) = A[G{x) sin 2 C + B(x) cos 2 <], 



(a; - a\f 



G(x) = exp I — 



2<7 2 



(Ti/2)^ 

(r 2 /2)^ 

(x-a 2 )ft! + (r 2 /2)ft! 



°y x ) - \ (r 2 /2)fe 



if x < ai\ 

if X > (22, 



(6) 



(7) 



(8) 



A, £, oi, cr, a2, Ti, f3i, T2, P2 are free fit parameters. The 
B(x) term is added to the Gaussian function to describe 
the asymmetric power-law tails of the 6m distributions. 

When used in data, the resolution a in the Gaussian 
term of Eq.©, is modified to take into account a possible 
difference between data and simulation: 



'MC 
'MC 



Act 2 if Act < 0; 
Act 2 if Act > 0. 



(9) 
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FIG. 9: The distribution of the difference (8m) between measured and true KsK ± ir^ 1 mass for simulated events (a) for the 
two-photon process e + e~ — ¥ e + e~rj c , and (b) for the ISR process e + e~ — > J/ip'). The curves correspond to fits defined in the 
text. 



The parameter ctmc is found to be 7.8 MeV/c 2 for the 
J/ip and 7.6 MeV/c 2 for the r/ c . The parameter Act is 
determined from the fit to the measured K$K it^ mass 
spectra. 

The background distribution is described by a second- 
order polynomial. Both spectra shown in Fig. [5] are fit 
simultaneously with 14 free parameters: the J/ip peak 
position (mj/^,), the difference between the J/ip and rj c 
mass values (Am), the rj c width (T), the numbers of r) c 
and J/ip events, Ac, and the background parameters for 
both spectra. The fitted curves are shown in Fig. [T0l 
For the full mass range, 2.7-3.3 GeV/c 2 , the \ 2 values 
corresponding to the r\ c and J/ip distributions (solid and 
dashed in Fig. [8]) are 230 and 198, respectively, for a total 
number of degrees of freedom 2 x 240 — 14. The resulting 
parameter values are as follows: 



Am 
F 

N no 
Act 



114.7 ±0.4 MeV/c 2 , 
31.7 ± 1.2 McV, 
14450 ± 320, 
-0.4 ±2.5 MeV/c 2 , 
3095.8 ±0.3 MeV/c 2 . 



(10) 



The mass resolution for the J/ip in data is found to be 
consistent with the prediction of MC simulation. The fit- 
ted value of the J/ip mass is shifted by — 1.1±0.3 MeV/c 2 
from the nominal J/ip mass value 0] ■ It is assumed that 
this mass scale shift does not affect the mass difference 
Am. Since the momentum distributions for J/ip and rj c 
events are different, and the MC simulation of the de- 
tector response is not perfect, we test this assumption as 
follows. The no-tag event sample was divided into three 
subsamplcs with approximately equal statistics but with 
different laboratory z momentum (p z ) of r\ c candidates. 
The average p z values in the subsamplcs are 3.2, —0.4, 



and —1.3 GeV/c, while the J/ip momentum is peaked 
at p z = —2.34 GeV/c. The fitted values of the Am pa- 
rameter for the three subsamplcs arc found to be shifted 
relative to the nominal fit value by 0.5 ± 0.6, —0.6 ± 0.6, 
and 0.2 ± 0.6 MeV/c 2 , respectively. We do not observe 
any significant dependence of the Am parameter on the 
r/c momentum direction and absolute value. Neverthe- 
less, the shift value at the maximum difference between 
the r/ c and J/ip momenta, 0.5 MeV/c 2 , is taken as an 
estimate of the Am systematic uncertainty due to the 
difference of the J/ip and r\ c momentum distributions. 

To estimate the uncertainty of the fit parameters due 
to the assumed background shape, the second-order poly- 
nomial describing background in Fig. llOf a) is replaced by 
an exponential function, and the changes in the param- 
eter values are considered to be measures of their asso- 
ciated systematic uncertainties. This yields AT = 0.8 
MeV and AiV^ = 400. 

The MC simulation uses a phase space distribution 
for ?y c — > KgK ± ^ decay. This can lead to distor- 
tion of the resolution function and a systematic change 
in the detection efficiency determined from simulation. 
In order to address this issue, a study of the Dalitz 
plot distribution for r/ c — > KsK^n^ decay was per- 
formed. The Dalitz plots for data events from the r\ c 
signal (2.94 < M KKn < 3.02 GeV/c 2 ) and sideband 
(2.90 < M KK „ < 2.94 GeV/c 2 and 3.02 < M kRit < 3.06 
GeV/c 2 ) regions are shown in Fig. 111! and their projec- 
tions are shown in Fig. [T2"l 

Figure [TBI shows the distributions of cos8k, where 9 k 
is the charged-kaon polar angle in the KgK^^ rest 
frame, for events from the r\ c region after background 
subtraction, and for events from the sidebands with the 
r/c contribution subtracted. The cos 9k distribution for 
r/c events closely resembles the distribution in the signal 



11 



^ 1000 
o 
> 

2 800 
in 



a 

> 



600 



400 



200 



: (a) 




1 

,i 
|,l 


l 


Jut 




'/ \ 


- 
- 




V 








1 


i 


2.8 




3 


3.2 



M KKjt (GeV/c1 




M^CGeV/O 



FIG. 10: The ^#±71^ invariant mass distribution and fitted curve for no-tag data events that (a) fail Eq. ([4j| , and (b) satisfy 
Eq.Q for ISR event selection. 
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FIG. 11: The Dalitz plots for data events from (a) the r\ c signal region 2.94 < M K ^ n < 3.02 GeV/c 2 , and (b) the combined 
sideband regions, 2.90 < M KRn < 2.94 GeV/c 2 and 3.02 < M kRw < 3.06 GeV/c 2 . 



MC simulation shown by the dashed histogram in FigfT^l 
Since the MC simulation uses a phase space decay 

model, this suggests that in the KsK ± tt t mode the r/ c 

decays predominantly via the scalar Xg(1430) meson, 

i.e., rj c -> KZ(U30)K. 

It should be noted that a significant part of the non-?7 c 

background comes from 77 — > KsK^ir^ , and so can in- 



terfere with the signal process 77 



Vc 



Ks^TT 



From the Dalitz plot, the background process seems 
to proceed mainly via the intermediate K*(892)K and 
K*(U30)K states, where the A* (1430) may be either 
the scalar (7^(1430)) or tensor (#2(1430)) state. In- 
terference between the 1 = and 1 = 1 background 
amplitudes appears to lead to the suppression of neutral 



K* (892) production (observed in Fig. H"27 a) with respect 
to Fig. H2T b)). The r\ c signal distribution in Fig. [13] ac- 
tually represents the cos 9 k dependence of the detector 
acceptance. The cos 6k distribution for non-?7 c events 
corrected for this acceptance will be peaked near ±1. 
This means that the kaon in continuum is dominantly 
produced with large orbital momentum. Therefore, in- 
terference between the ^-wave r\ c decay amplitude and 
the 77 — > KsK^-rT nonresonant amplitude is expected 
to be small. 



To estimate possible shifts of the rj c parameter values 
due to interference, an interference term is introduced 
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FIG. 12: The projections of the Dalitz plots of Fig. 1111 the solid and dashed histograms correspond to the r) c signal and 
sideband M K g^ regions, respectively. 




cos 9 



FIG. 13: The cosQk distribution for events from the r\ c signal 
region after background subtraction (points with error bars), 
and for events from the sidebands after subtraction of the r\ c 
contribution (solid histogram). The dashed histogram repre- 
sents the signal MC simulation. 



GeV 2 /c 4 , where in both signal and background the quasi- 
two-body final states, K*(U30)K and K%(U30)K, dom- 
inate. Therefore, no significant phase shift is expected 
between the two amplitudes, and so parameter A is cho- 
sen to be real. From the fit, A = 0.03 ± 0.01 and there 
are insignificant changes in the values of parameters T 
and N Vc . The value of Am changes by 1.5 MeV/c 2 , and 
so this is considered to provide an estimate of the sys- 
tematic uncertainty due to possible interference between 
the r) c and nonresonant amplitudes. 

To take into account the difference between data and 
simulation of the r\ c — > KsK^ir^ decay dynamics, the 
Dalitz plot is divided into 26 cells as shown in Fig. [TTT a). 
For each cell the KsK ± tt t mass spectrum is fit using 
the fitting function described above, and the number of 
rj c events is determined. This experimental Dalitz plot 
distribution corrected for detection efficiency is used to 
reweight the signal simulation. The reweighting changes 
the resolution function and the fit parameter values in- 
significantly. 

Thus, from the fit to the KsK ± tt t mass spectrum 
for no-tag events the following r\ c parameter values are 
obtained: 



into the fitting function through the following: 



M -M -iT/2 +A y P 2 (M ) ' 

where the Brcit-Wigner function describes the rj c ampli- 
tude, P2(M) is a second-order polynomial describing the 
mass dependence of the nonresonant intensity, and A is 
the value of this amplitude at the r\ c mass. The Pz(M) 
coefficients are chosen to be equal to the coefficients of 
the second order polynomial describing the nonresonant 
background. A comparison of the Dalitz plot distribu- 
tions for r\ c and non-?7 c events indicates that the maximal 
interference should occur in the vicinity of M 2 (Ktt) « 2 



Am = 114.7 ± 0.4 ± 1.6 MeV/c 2 , 

T = 31.7 ± 1.2 ±0.8 MeV, 
N Vc = 14450 ± 320 ±400. (12) 

When the nominal value of the J/ip mass Q is used, 
the T) c mass becomes 2982.2 ± 0.4 ± 1.6 MeV/c 2 . The 
results for the mass and width are in agreement with 
the previous BABAR measurement obtained using 88 fb" 1 
data [ll|: rn Ve = 2982.5 ± 1.1 ± 0.9 MeV/c 2 and T = 
34.3±2.3±0.9 MeV. The systematic uncertainty reported 
here is greater than that reported in the previous analy- 
sis as we are now allowing for the possibility that there 
exists a J p = 0~ continuum KsK ± TT zfl amplitude which 
interferes with the rj c amplitude. 
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FIG. 14: The KsK ± tt^ 1 mass spectrum for single-tag data 
events with 2 < Q 2 < 50 GeV 2 . The solid curve is the fit 
result. The dashed curve represents non-peaking background. 



V. FITTING THE KsK*** MASS SPECTRUM 
FOR SINGLE- TAG EVENTS 

The KsK ± tt t mass spectrum for single-tag events 
from data with 2 < Q 2 < 50 GeV 2 is shown in Fig. [TJ] 
For Q 2 > 50 GeV 2 we do not see evidence of an r\ c signal 
over background. To determine the number of rj c events, 
a binned likelihood fit is performed to the spectrum using 
a function consisting of a sum of r] c , J/ip, and background 
distributions. The bin size used in the fit is chosen to be 
2.5 MeV/c 2 . The mass resolution line shape is described 
by the function of Eq. ([6]) with parameters determined 
from the signal simulation reweighted to reproduce the 
Q 2 dependence observed in data. This resolution func- 
tion and its convolution with a Breit-Wigner function are 
used to describe the J/ip and rj c line shapes, respectively. 
The background distribution is described by either a sec- 
ond order polynomial or an exponential function. The 
fit result is shown in Fig. [T3]for a quadratic background. 
The fitted r] c parameter values, Am = 111.2±2.0 MeV/c 2 
and r = 31.9 ± 4.3 MeV, are in agreement with the re- 
sults obtained for no-tag events, and the total r] c signal 
is 530 ± 41 events. The difference in signal yield for the 
two background hypotheses is 17 events. 

A fitting procedure similar to that described above is 
applied in each of the eleven Q 2 intervals indicated in 
Table HI The parameters of the mass resolution function 
are taken from the fit to the mass spectrum for simu- 
lated events in the corresponding Q 2 interval. The r\ c 
and J/tp masses are fixed at the values obtained from 
the fit to the spectrum of Fig. [TH while the r\ c width is 
taken from the fit to the no-tag data. The free parame- 
ters in the fit are the numbers of r\ c and J/ip events, and 
two or three additional parameters depending upon the 



description of the background shape (quadratic or expo- 
nential). The KsK^ir^ mass spectra and fitted curves 
(quadratic background) for three representative Q 2 inter- 
vals are shown in Fig. 1151 The number of rj c events ob- 
tained from the fit using a quadratic background is listed 
for each Q 2 interval in Table HI The difference between 
the fit results for the two background hypotheses is used 
as an estimate of the systematic uncertainty associated 
with the assumed background form. 



VI. PEAKING BACKGROUND ESTIMATION 
AND SUBTRACTION 

Background containing true ?7 c 's might arise from e + e~ 
annihilation processes and two-photon processes with 
higher multiplicity final states. The processes with a J/tp 
in the final state considered in previous sections are also 
sources of peaking background because of the relatively 
large branching fraction for the decay J/ijj — > r\ c ^ 0). 

For no-tag events the most discriminating variable be- 
tween signal and background is the r\ c candidate trans- 
verse momentum (p* ± )- In particular it is expected that 
background from e + e~~ annihilation increases rapidly 
with transverse momentum. Figure [TBI shows the p*j_ dis- 
tribution for no-tag data events containing an r\ c . In 
each p* ± interval the number of rj c events is determined 
from the fit to the KsK^n^ mass spectrum. It is seen 
that the data distribution is in good agreement with sig- 
nal MC simulation. A conservative upper limit on the 
level of e + e~ annihilation background is obtained by fit- 
ting a sum of the MC signal distribution and a constant 
background to the data p* ± distribution. The number of 
background events with p* ± < 0.25 GeV/c is found to be 
110 ± 150. 

The two-photon background from the process e + e~ — > 
e + e~?/ c 7r° is studied using a special selection. From the 
sample of events satisfying preliminary selection crite- 
ria, events with two or more extra photons are selected 
with the energy of each photon required to be greater 
than 50 MeV. Two photons with total energy greater 
than 0.2 GeV, and invariant mass in the range 0.07-0.20 
GeV/c 2 form a ir° candidate. The 7r° candidate is com- 
bined with an ?y c candidate, and it is required that the 
transverse momentum of the r/ c 7r° system be less than 
0.25 GeV/c. The resulting tt° candidate mass spectrum 
is shown in Fig. [171 and the KsK^-k^ mass spectrum 
for events with 0.11 < M 77 < 0.15 GeV/c 2 is shown in 
Fig. [TSJ The spectrum is fitted by a sum of signal (rj c 
and J/ip) and background functions. The fitted num- 
ber of T] c events is found to be 60 ± 40. Simulation 
of the e + e~ — > e + e~rj c Tr° process shows that the ra- 
tio of the numbers of events selected using the standard 
and special criteria is about 2.5, so that the estimated 
background from e + e~ — > e + e~rj c TT° in the event sam- 
ple with the standard selection is 150 ± 100. A similar 
approach is used to estimate a possible background from 
the e + e~ — > e + e~rj c rj process. No rj signal is observed in 
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FIG. 15: The KsK^if^ mass spectra for single-tag data events from three representative Q 2 intervals. In each figure the solid 
curve represents the fit result, and the dashed curve indicates non-peaking background. 



TABLE I: The Q 2 interval, number of events with r/ c obtained from the fit (A^ c ), number of background events from J/tp — > r\ c ^ 
decay (iVbkg), efficiency correction (5totai), and number of signal events corrected for data/MC difference and resolution effects 
(iVcTrr° ldcd )- The quoted errors on N Vc and N^° ldcd are statistical and systematic. 



Q 2 interval (GeV 2 ) 


N Vc 


A'bkg 


5total(%) 


TV ru ii folded 
* * corr 




no-tag 


14450 ± 320 ± 400 


730 ± 240 


-2.6 


14090 ± 330 ± 480 


14.5 


2-3 


41.0 ± 8.6 ± 1.3 


0.7 ±0.4 


-0.1 


39.9 ±9.0 ±1.4 


2.0 


3-4 


56.2 ±10.5 ±4.0 


0.6 ±0.4 


0.0 


55.3 ±10.9 ±4.2 


4.9 


4-5 


65.0 ±10.9 ±1.1 


0.1 ±0.4 


-0.1 


64.8 ± 11.5 ± 1.2 


9.1 


5-6 


52.6 ±9.6 ±0.6 


0.5 ±0.4 


-0.4 


51.8 ± 10.3 ±0.8 


12.1 


6-8 


90.9 ± 12.2 ±4.6 


1.3 ±0.8 


-0.4 


90.3 ± 12.8 ±4.9 


14.0 


8-10 


60.9 ±10.9 ±2.8 


0.9 ±0.6 


-0.8 


61.3 ± 11.7 ±3.1 


17.9 


10-12 


34.8 ± 7.3 ± 1.8 


1.0 ±0.6 


-1.0 


33.5 ±7.9 ±2.1 


21.4 


12-15 


42.3 ±8.7 ±2.1 


1.9 ±0.8 


-1.3 


41.2 ±9.4 ±2.4 


23.0 


15-20 


45.5 ±7.9 ±1.0 


2.4 ± 1.0 


-1.0 


44.3 ±8.5 ±1.5 


23.8 


20-30 


23.7 ±6.6 ±0.6 


1.6 ±0.7 


-1.0 


22.5 ±6.9 ±1.0 


24.7 


30-50 


10.8 ±4.5 ±0.1 


0.9 ±0.5 


-1.3 


10.3 ±4.8 ±0.5 


24.5 




0.1 0.2 0.3 0.4 0.5 0.1 0.15 0.2 

P± (GeV/c) M (GeV/c 2 ) 



FIG. 16: The distribution of the r/ c candidate transverse mo- FIG. 17: The tt° candidate mass spectrum for the selected 
mentum for no-tag data events containing an r\ c (points with e + e~ — \ e + e~ K,sK ± tt^ 1 tt data events in no-tag mode, 
error bars), and for simulated signal events (histogram). 




FIG. 18: The KsK if* mass spectrum for e + e — > 
e + e~ KsK ± tt^ 7r° candidate events in no-tag mode with 
0.11 < M 77 < 0.15 GeV/c 2 . The solid curve is the fit re- 
sult. The dashed curve represents non-peaking background. 



FIG. 19: The detection efficiency as a function of momentum 
transfer squared for events with a tagged electron (squares), 
a tagged positron (triangles), and their sum (circles). 



the two-photon mass spectrum, nor is there an r/ c signal 
in the KgK^n^ mass spectrum. The e + e~ — > e + e~r\ c r] 
background is therefore considered negligible. 

Background from both sources, e + e~ annihilation and 
two-photon processes, does not exceed 490 events (90% 
CL), i.e., 3.5% of the total number of r\ c events. This 
number is considered to provide an estimate of the sys- 
tematic uncertainty due to possible e + e~ annihilation 
and two-photon background. 

The total number of e + e~ — > J/tpJ, J/i/> —> 
KsK^tt^ events found in the no-tag event sample is 
3170 ± 100. From simulation the ratio of the detec- 
tion efficiencies e(e + e~ — > J/V'7, J/ip — > ?7c7j Vc — >• 
KsK ± n T )/e(e + e~ -4 J/tpf, J/i> -> KsK^iC?) is found 
to be 1.18 ± 0.01. Taking into account the ratio of the 
branching fractions 



b = 



B(J/4> -> KKtt) 



0.20 ±0.07 (13) 



the corresponding background contribution to the rj c 
peak from the ISR process is found to be 730±240 events. 

For single-tag events, the background from e + e~ an- 
nihilation can be estimated using events with the wrong 
sign of the e^r\ c momentum z-component. The KsK^tt^ 
mass spectrum for the wrong-sign events is shown in 
Fig. [7] together with the spectrum for right-sign events. 
Assuming that the numbers of background events for the 
wrong- and right-sign data samples are approximately 
the same, a fit to the mass spectrum for wrong-sign events 
yields 1.4 ± 3.0 e + e~ annihilation events peaking at the 
rjc mass. 

To estimate the two-photon background in the single- 
tag event sample, the result obtained for no-tag events is 



used. Assuming that the signal and background Q 2 de- 
pendences are approximately the same, it is estimated 
that the number of two-photon background events is 
5.7 ± 3.0. The total background from e + e~ annihila- 
tion and two-photon processes does not exceed 13 events 
(90% CL), or 2.5% of the total number of r/ c events in 
the single-tag event sample. This number is considered to 
provide an estimate of the systematic uncertainty due to 
possible e + e~ annihilation and two-photon background. 

The background from the process e + e~ — > e + e~J/%jj 
(Fig. Hb)), J/tp -> tj c j, i lc -> KsK^it* is esti- 
mated from the measured Q 2 distribution (Nj/^,^) for 
the e + e~ — > e + e~J/tp, J/?/' — > KsK^tt^ events as 



ibN 



where b is the ratio of the branching frac- 



tions defined in Eq. (fT3]) , and m is the ratio of the detec- 
tion efficiencies for the J/ip —> ?/ c 7, r\ c — > KsK^tt 11 and 
J/tp —> KsK ±r ir^ decay modes. The coefficient m varies 
from 0.7 to 0.5 in the Q 2 range of interest. The estimated 
background contributions resulting from J/tp — > rjcj de- 
cay are listed in Table Q] The fraction of background 
in the e^rjc data sample changes from about 1.0% for 
Q 2 < 10 GeV 2 to about 5% at Q 2 w 30 GeV 2 . 



VII. DETECTION EFFICIENCY 

The detection efficiency is determined from MC simu- 
lation as the ratio of the true Q 2 distributions computed 
after and before applying the selection criteria. The Q 2 
dependence of the detection efficiency is shown in Fig.H^l 
The detector acceptance limits the detection efficiency at 
small Q 2 . The cross section is measured in the region 
Q 2 > 2 GeV 2 where the detection efficiency is greater 
than 2%. The asymmetry of the e + e~ collisions at PEP- 
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FIG. 20: The Q 2 dependence of the relative difference be- 
tween detection efficiencies determined from MC simulation 
with and without Dalitz-plot reweighting. 
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FIG. 22: The correction to the MC-estimated kaon identifi- 
cation efficiency as a function of Q 2 . 
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FIG. 21: The Q 2 dependence of the data-MC double ratio for 
events with and without the polar angle restriction 0.387 < 
6 < 2.400 rad for all four charged-particle tracks from the r) c 
decay. 



II leads to different efficiencies for events with electron 
and positron tags. The Q 2 range from 2 to 6 GeV 2 is 
measured only with the positron tag. For no-tag events 
the detection efficiency is 0.1446 ± 0.0023. The efficiency 
is calculated using simulated events reweighted accord- 
ing to the Dalitz plot distribution observed in data. For 
no-tag events the relative difference between detection ef- 
ficiencies calculated with and without weighting is found 
to be —(1.1 ± 1.6)%. The quoted error is determined by 
the statistical errors of the measured Dalitz plot distribu- 
tion. The corresponding relative difference for single-tag 
events is shown in Fig.J^Hlas a function of Q 2 . The detec- 
tion efficiency has only a very weak dependence on the 
dynamics of rj c decay. 

Possible sources of systematic uncertainty due to differ- 
ences between data and simulated detector response are 
now considered. For no-tag events the MC simulation 



FIG. 23: The correction to the MC-estimated electron iden- 
tification efficiency as a function of Q 2 . 



predicts a significant loss of signal events, (5.3 ± 0.1)%, 
due to background filters used in event reconstruction. 
The filter inefficiency can be measured in data using a 
special sample of prescaled events that does not pass 
the background filters. The filter inefficiency obtained 
in data is (7.5 ±1.2)%. The difference 5 = -(2.2 ±1.2)% 
is used to correct the number of signal events. For single- 
tag events the presence of the additional electron leads to 
a significantly smaller filter inefficiency. The simulation 
predicts a filter inefficiency of (0.57 ± 0.02)%, which is 
about ten times smaller than for no-tag events. We con- 
clude that this source of systematic uncertainty is negli- 
gible for single-tag events. 

To study the possible systematics for no-tag events due 
to selection criteria, the mass window for the K$ can- 
didate is increased from 487.5-507.5 MeV/c 2 to 475.0- 
520.0 MeV/c 2 , the limit on the rj c candidate transverse 
momentum is changed from 0.25 to 0.50 GeV/c, and the 
polar angle restriction is set to 0.387 < 6 < 2.400 rad for 
all four charged-particle tracks in an event. The last mod- 
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TABLE II: The main sources of systematic uncertainty as- 
sociated with the detection efficiency and the total efficiency 
systematic error for no-tag and single-tag events. 



Source 


No-tag, % 


Single-tag, % 


trigger, filters 


1.2 




rj c selection 


5.9 


5.7 


track reconstruction 


1.4 


1.5 


K identification 


0.4 


0.5 


e identification 




0.5 


total 


6.2 


5.9 



ification rejects about 30% of signal events. The double 
ratio 



(N nev ,/N) 



data 



(iV new /iV)MC 



(14) 



is calculated, where N nevf and N are the numbers of sig- 
nal events with the new and standard selection criteria, 
and is found to be close to unity for the definition of the 
Ks mass window (0.993 ± 0.005), and for the condition 
onp* (1.002 ± 0.009). A significant deviation from unity 
(5.9 ± 1.8)% is observed for the polar angle restriction. 
This deviation is taken as an estimate of the systematic 
uncertainties due to imperfect simulation of the r\ c selec- 
tion criteria close to the limits of the fiducial tracking 
region. 

The systematic uncertainty due to a possible difference 
between the data and simulation in the charged-particle 
track reconstruction for pions and kaons is estimated to 
be about 0.35% per track, so that the total uncertainty 
is estimated to be 1.4%. The data-MC simulation dif- 
ference in the kaon identification is estimated using the 
identification efficiencies measured for kaons from the 
D* + D°7r+, D° -> ir + K~ decay sequence. The ef- 
ficiency correction is found to be —0.4%. The systematic 
uncertainty associated with this correction is taken to 
be equal to the value of the correction, 0.4%. The to- 
tal efficiency correction for no-tag events is —2.6%, and 
the total systematic error associated with the efficiency 
is about 6.2%. 

The polar angle restriction described above is also 
tested for single-tag events, and Fig. [21] shows the Q 2 
dependence of the data-MC double ratio. No signifi- 
cant Q 2 dependence is observed, and the average value is 
1.057 ± 0.032, which is very close to the value for no-tag 
events. For single-tag events the systematic uncertainty 
due to the rj c selection criteria is estimated to be 5.7%. 
The efficiency correction (<5i) for kaon identification for 
single-tag events is shown in Fig. [22] function of Q , 
and this results in an associated systematic uncertainty 
of 0.5%. 

The data-MC simulation difference in electron identi- 
fication is estimated using the identification efficiencies 
measured for electrons from radiative Bhabha events. 
The efficiency correction ($2) is shown as a function of 
Q 2 in Fig. [53J The associated systematic uncertainty 



does not exceed 0.5%. The systematic uncertanty due to 
data-MC simulation difference in the electron track re- 
construction is about 0.1%. To estimate the effect of the 
requirement —0.02 < r < 0.03 (see Fig. [5]), events with 
0.03 < r < 0.06 are studied. The data-MC simulation 
double ratio defined by Eq. fT^)) is found to be consistent 
with unity, (0.99 ± 0.02), so that the simulation repro- 
duces the shape of the r distribution very well. 

The main sources of systematic uncertainty associ- 
ated with detection efficiency are summarized in Table [TT] 
for the no-tag and single-tag samples. The values of 
the detection efficiency and total efficiency correction 
<5totai,i ~ Si t i + 62,1 for different Q 2 intervals are listed 
in Table [U The data distribution is corrected as follows: 



Ni/(1 



•'total, i 



(15) 



where Ni = N n 
in the ith Q 2 interval. 



-^bkg.i is the number of signal events 



VIII. CROSS SECTION AND FORM FACTOR 



The Born differential cross section for 
f e~?7 c ,77 c — > KKtt is 



dcr 

dQ 2 



(dN/dQ 2 ) 



2 "\ unfolded 
corr 



eRL 



B{r} c -> KKtt) 



B(j] c -> K s K±TrT)B(K s 



(16) 



where (dA^/d<5 2 )™ r f ° ldod is the mass spectrum corrected 
for data-MC simulation difference and unfolded for de- 
tector resolution effects, explained below, L is the total 
integrated luminosity, e is the Q 2 -dependcnt detection 
efficiency, R is a radiative correction factor accounting 
for distortion of the Q 2 spectrum due to the emission of 
photons from the initial-state particles, and for vacuum 
polarization effects. The ratio B(r/ C — s- KKtt) /B(r) c — ► 
KsK^nT) is expected to be 3 from isospin relations. 

The radiative correction factor is determined using 
simulation at the generator level (with no detector sim- 
ulation). The Q 2 spectrum is generated using only the 
pure Born amplitude for the e + e~ — > e + e~rj c process, 
and then using a model with radiative corrections in- 
cluded. The Q 2 dependence of the radiative correction 
factor, evaluated as the ratio of the second spectrum to 
the first, is shown in Fig.(2H and is fitted using the func- 
tion a/ (1+6Q 7 ). The accuracy of the radiative correction 
calculation is estimated to be 1% @. It should be noted 
that the value of R depends on the requirement on the 
extra photon energy. The Q 2 dependence obtained cor- 
responds to the condition r = 2E*/t/s < 0.1 imposed in 
the simulation. 

The corrected and unfolded Q 2 distribution 
(diV/d(3 2 )cor f ° ldcd is obtained from the measured 
distribution by dividing by the efficiency correction 
factor (see Eq.(Qj)])) and unfolding the effect of Q 2 
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FIG. 24: The Q 2 dependence of the radiative correction fac- 
tor. 



resolution. Using MC simulation, a migration matrix 
H is obtained, which represents the probability that 
an event with true Q 2 in interval j be reconstructed in 
interval i: 



dN 
dQ 2 



dN 
dQ 5 



(17) 



In the case of extra photon emission, Q 2 ruo i s calculated 
as —{p — p 1 — k) 2 , where k is the photon four-momentum; 
e and R in Eq. (fl"6"|) are functions of Q 2 ruc - As the cho- 
sen Q 2 interval width significantly exceeds the resolution 
for all Q 2 , the migration matrix is nearly diagonal, with 
values of diagonal elements ~ 0.95, and of the next-to- 
diagonal ~ 0.02. The true Q 2 distribution is obtained 
by applying the inverse of the migration matrix to the 
measured distribution. The procedure changes the shape 
of the Q 2 distribution insignificantly, but increases the 
errors (by ~10%) and their correlations. The corrected 
Q 2 spectrum (V c u n r f° ldod ) is listed in Table U 

The values of the differential cross sections obtained 
are listed in Table IIII1 where the first error is statis- 
tical and the second systematic. The latter includes 
only Q 2 -dcpcndent errors, namely, the systematic uncer- 
tainty in the number of signal events and the statistical 
errors on the efficiency correction and MC simulation. 
The Q 2 -independent systematic error is 6.6%; this re- 
sults from the systematic uncertainties on the detection 
efficiency (5.9%), background subtraction (2.5%), the ra- 
diative correction factor (1%), and the error on the in- 
tegrated luminosity (1%) combined in quadrature. The 
MC simulation for single-tag events is performed, and 
the detection efficiency is determined, with the restric- 
tion that the momentum transfer to the untagged elec- 
tron be greater than —1 GeV 2 , so that the cross section 
is measured for the restricted range \q 2 \ < 1 GeV 2 . The 
measured differential cross section is shown in Fig. [231 

Because of the strong nonlinear dependence of the 
cross section on Q 2 , the value of Q 2 corresponding to the 



TABLE III: The Q 2 interval and the weighted average Q 
value (Q 2 ), the e + e~ — > e + e~r] c cross section multiplied by 
B(rj c — > KKn) [da/dQ 2 (Q 2 )], and the normalized 77* — > r\ c 
transition form factor (\F(Q 2 ) /F(0)\). The statistical and 
systematic errors are quoted separately for the cross section, 
but are combined in quadrature for the form factor. Only Q 2 - 
dependent systematic errors are quoted; the Q 2 -independent 
error is 6.6% for the cross section and 4.3% for the form factor. 



Q 2 interval 


Q 2 


da/dQ 2 (Q 2 ) 
(fb/GeV 2 ) 


\F(Q 2 )/F(0)\ 


(GeV 2 ) 


(GeV 2 ) 




2-3 


2.49 


18.7 ±4.2 ±0.8 


0.740 ± 0.085 


3-4 


3.49 


10.6 ±2.1 ±0.8 


0.680 ± 0.073 


4-5 


4.49 


6.62 ±1.18 ±0.19 


0.629 ± 0.057 


5-6 


5.49 


4.00 ±0.80 ±0.10 


0.555 ± 0.056 


6-8 


6.96 


3.00 ±0.43 ±0.17 


0.563 ± 0.043 


8-10 


8.97 


1.58 ±0.30 ±0.08 


0.490 ± 0.049 


10-12 


10.97 


0.72 ±0.17 ±0.05 


0.385 ± 0.048 


12-15 


13.44 


0.55 ±0.13 ±0.03 


0.395 ± 0.047 


15-20 


17.35 


0.34 ±0.07 ±0.01 


0.385 ± 0.038 


20-30 


24.53 


0.084 ± 0.026 ± 0.004 


0.261 ±0.041 


30-50 


38.68 


0.019 ± 0.009 ± 0.001 


0.204 ± 0.049 



measured cross section differs slightly from the center of 
the Q 2 interval. The measured cross section is described 
by a smooth function, which is then used to reweight the 
simulated Q 2 distribution and calculate the weighted av- 
erage value (Q 2 ) for each Q 2 interval. The values of Q 2 
obtained are listed in Table IIIII 

The no-tag event sample is used to obtain the total 



cross section for the reaction e H 
KKtt: 



~Vc, Vc 



Nr., 



Ctotal 



tL 



B(rj c -> KKtt) 



B( Vc -> K s K±ttT)B(K s 



(18) 



In the no-tag mode the radiative correction is expected 
to be less than 1% [l2| . and so the associated system- 
atic uncertainty is assigned conservatively as 1%. Taking 
V corr and e from Table HI the value obtained from Eq. (fT8]) 
is 



ototai = 0.900 ± 0.021 ± 0.074 pb. 



(19) 



The quoted errors are statistical and systematic, respec- 
tively. The latter is 8.1% and includes the systematic un- 
certainty in the number of signal events (3.3%), the sta- 
tistical and systematic errors on the detection efficiency 
(1.6% and 6.2%, respectively), on the background sub- 
traction (3.5%), on the radiative correction (1%), and 
the error on the integrated luminosity (1%). Using MC 
simulation, the calculated total cross section correspond- 
ing to r(?y c — > 77)£>(?y c — > KKtt) = 1 keV is found to be 
2.402 pb, and hence from Eq.(fT9"|) the value 



T(r) c -4 "f"f)B{rjc 



KKtt) = 
0.374 ± 0.009 ± 0.031 keV (20) 
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FIG. 25: The e + e" -> e + e"7? c differential cross section mul- 
tiplied by the r] c — > KKn branching fraction; the statistical 
and Q 2 -dependent systematic errors of Table Mil have been 
combined in quadrature. 



FIG. 26: The 77* — ¥ rj c transition form factor normalized to 
F(0) (points with error bars). The solid curve shows the fit 
to Eq. (|22p . The dotted curve shows the leading order pQCD 
prediction from Ref. §]■ 



is obtained. This result agrees with the Particle Data 
Group value 0.44 ± 0.05 keV Q, and also with the recent 
CLEO measurement 0.407 ± 0.022 ± 0.028 keV 0. 

To extract the transition form factor, the measured 
and calculated Q 2 spectra are compared. The simulation 
for single-tag events uses a constant form factor value, so 
that the measured normalized form factor is determined 
from 



|F 2 (Q 2 )/F 2 (0)| 



(d7V/dQ 2 )^ lc . tag 

Ardata 
no-tag 

- _MC 
no-tag^ total 



£ s i„ g ic-ta g (g 2 )(da/dQ 2 ) s R f c 



(21) 



single-tag 



The normalized form factor is proportional to the ra- 
tio of the number of single-tag events to the number of 
no-tag events. It is expected that part of the systematic 
uncertainty, in particular that associated with the detec- 
tion efficiency, cancels in this ratio. However, the single- 
tag data sample is not large enough for a detailed study 
of the Q 2 -dependence of the observed data-MC simula- 
tion difference in detector response. Therefore, this un- 
certainty is estimated conservatively to be equal to the 
corresponding systematic uncertainty for the single-tag 
events (6%). The values of the form factor obtained are 
listed in Table IIII1 and shown in Fig. [211 with the sta- 
tistical and Q 2 -dependcnt systematic errors combined. 
The <5 2 -independent systematic error on the form factor 
is 4.3%; this value combines in quadrature the system- 
atic uncertainty on detection efficiency, the uncertainty 
on the number of no-tag events, the statistical error on 
the detection efficiency for no-tag events, the uncertain- 
ties associated with the background subtraction, and the 
uncertainty on the radiative correction. 



The form factor data of Fig. 
the monopole form 

\F(Q 2 )/F(0 



I are well described by 



(22) 



1 + Q 2 /A' 



as shown by the solid curve. The corresponding fitted 
value of the pole parameter A is 



A = 8.5 ±0.6 ±0.7 GeV 2 , 



(23) 



where the second quoted error is due to the 4.3% Q 2 - 
independent systematic error on the measurements. This 
value of the pole parameter is in reasonable agreement 
with that expected from vector dominance, namely A = 
m 'j/ip = 9-6 GeV 2 , and in good agreement with the lat- 
tice QCD calculation, A = 8.4 ± 0.4 GeV 2 @. The dot- 
ted curve in Fig. |2"61 shows the result of the leading-order 
pQCD calculation of Ref. Q. The data lie systematically 
below this prediction, but within the theoretical uncer- 
tainty quoted in Ref. [||. 



IX. SUMMARY 

The reaction e + e~ — > e + e~r] c , with T) c — > KsK^ir^ , 
has been studied in the no-tag and single-tag modes. We 
measure the following values for the r\ c mass and width: 

m, c = 2982.2 ± 0.4 ± 1.6 MeV/c 2 , 

T = 31.7 ± 1.2 ±0.8 MeV. (24) 

These results agree with earlier BABAR measurements [Tl| 
and supersede them. 
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We have also measured the total cross section 

a(e + e~ — > e + e~rj c )B(i] c — > KKir) and differential cross 

section (da/dQ 2 )B(r] c —> KKn). From these data we 
determine the value 

r (?7c -> li)B{r) c -> if -KV) = 

0.374 ± 0.009 ± 0.031 keV (25) 

and measure the normalized 77* — > transition form 
factor \F(Q 2 )/ F(0)\ for the momentum transfer range 
from 2 to 50 GcV 2 . The latter is well described by the 
simple monopole form of Eq. (f2"2"|) with A = 8.5 ± 0.6 ± 
0.7 GeV 2 in agreement with both the vector dominance 
expectation and the QCD prediction. 
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